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Since 1984, electric cell-substrate impedance sensing (ECIS) has been used to monitor cell behavior 
in tissue culture and has proven sensitive to cell morphological changes and cell motility. We have 
taken ECIS measurements on several cultures of non-cancerous (HOSE) and cancerous (SKOV) 
human ovarian surface epithelial cells. By analyzing the noise in real and imaginary electrical 
impedance, we demonstrate that it is possible to distinguish the two cell types purely from signatures 
of their electrical noise. Our measures include power-spectral exponents. Hurst and detrended 
fluctuation analysis, and estimates of correlation time; principal-component analysis combines all 
the measures. The noise from both cancerous and non-cancerous cultures shows correlations on 
many time scales, but these correlations are stronger for the non-cancerous cells. 

PACS numbers: 87.18.Ed, 87.80.Tq, 05.40.Ca 
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I. INTRODUCTION 

Electrical cell-substrate impedance sensing (ECIS) has 
been in use since 1984 [l| to monitor changes in cell cul- 
tures due to spreading or in response to chemical stim- 
uli, infection, or flow. Applications include studies of cell 
migration, barrier function, toxicology, angiogenesis, and 
apoptosis. Several papers have noted that impedance 
fluctuations are associated with cellular micromotion ^ . 
However, we are not aware of any previous work apply- 
ing statistical techniques to these fluctuations in order 
to distinguish two different cell types. Here, we demon- 
strate that measures of the electrical noise from cultures 
of cancerous and non-cancerous human ovarian surface 
epithelial cells distinguish them. We find that the noise 
in both cancerous and non-cancerous cultures shows cor- 
relations on many time scales, but by all measures, these 
correlations are weaker or of shorter range in the cancer- 
ous cultures. 



II. EXPERIMENTAL METHODS 

We used the ECIS system to collect micro-motion time- 
series data, the fluctuations in which are caused by the 
movements in a confluent layer of live cells. The system 
can be modeled as an RC circuit The cells 

are cultured on a small gold electrode (5 x 10~^cm^), 
which is connected in series to a 1-Megaohm resister, an 
AC signal generator operating at 1 volt and 4000 Hz, and 
flnally to a large gold counter-electrode (0.15 cm^). This 
network is connected in parallel to a lock-in amplifler, 
and the in-phase and out-of-phase voltages are collected 
once a second, from which we extract time series of resis- 
tance and capacitive reactance (Figure [1^) . In ECIS ex- 
periments, the fluctuations in complex impedance come 
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primarily from changes in intercellular gaps and in the 
narrow spaces between the cells and the small gold elec- 
trode 0) S @] ■ A current of about one microamp is driven 
through the sample, and the resulting voltage drop of a 
few millivolts across the cell layer has no physiological 
effect: this is a noninvasive, in vitro-technique. An ovar- 
ian cancer line (SK0V3) and a normal human ovarian 
surface epithehal (HOSE) cell hue (H0SE15) were pro- 
vided by Dr. Samuel Mok at Harvard Medical School. 
These cells were grown in M199 and MCDB 105 (1:1) 
(Sigma, St. Louis, MO) supplemented with 10% fe- 
tal calf serum (Sigma), 2mM L-glutamine, 100 units/ml 
penicillin, and 100 microgram/ml streptomycin under 
5% CO2, and a 37° C, high-humidity atmosphere. For 
ECIS micro-motion measurements, cells were taken from 
slightly sub-confluent cultures 48 hours after passage, 
and a mono-disperse cell suspension was prepared using 
standard tissue-culture techniques with trypsin/EDTA. 
These suspensions were equilibrated at incubator condi- 
tions before addition to the ECIS electrode wells. Con- 
fluent layers were formed 24 hours after inoculation, re- 
sulting in a density of lO^cell/cm^. 

Figure[T^ shows a representative 4096-second run (just 
over one hour) measuring the real part of impedance as 
a function of time; the example shows a HOSE culture, 
but to the eye, SKOV cultures do not appear very differ- 
ent. While the example shows increasing resistance with 
time, others show a decrease; at this time scale, there 
is no evidence for an overall trend. We collected, under 
similar conditions, 18 time series for HOSE cultures, of 
which 16 went for 8192 seconds and two for 4096 seconds. 
Each 8192-second run was split in two halves, so that ef- 
fectively we had thirty-four 4096-second runs; however, 
where appropriate in the analysis below, we discard the 
second halves of the longer runs in order to avoid inad- 
vertently introducing correlations. Similarly, for SKOV 
cultures we took data in eight 8192-second runs and ten 
4096-second runs, yielding effectively twenty-six 4096- 
second runs. We numerically differentiated the resistance 
and capacitance time series to obtain noise time series for 
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FIG. 1: Scheme of data extraction from noise, (a) Time series 
of resistance for one of the experimental runs. Taking the dis- 
crete time derivative and normalizing to zero mean and unit 
variance gives the noise, (b). The power spectrum of noise 
is shown in (c), using overlapping windows of 256 points in 
order to reduce scatter. Fits to the first hundred and last 
hundred frequencies estimate low- and high-frequency power- 
laws, / White noise would have appeared frequency- 
independent {a — 0). The Fourier transform of the power 
spectrum gives the autocorrelation, (d), which we fit to a 
shifted power- law decay and extract the measure l3o. As ex- 
plained in the text, subtle differences in the univariate noise 
distribution (e) (smoothed) discriminate between cancerous 
and non-cancerous micromotion. 

each, which we normalized to zero mean and unit vari- 
ance (Figure [TJd) . 

III. STATISTICAL MEASURES OF NOISE 

We seek information from the normaUzed noise series. 
The first question to pose is whether the noise can distin- 
guish cancerous from non-cancerous cultures, but more 
generally the measures we extract may be used to test 
models of cell micromotion. Broadly, such models may 
be characterized by short-term and long-term correlation, 
so we look at several measures for each. 

First, the power spectral density (Figure[TJ:) looks very 
much more like "pink noise" than "white noise;" that is. 



it shows signs of long-time correlations. A log-log plot 
of spectral density against frequency, /, suggests an in- 
tensity going as in the low- frequency limit. (We dis- 
cuss below the extent to which a true white-noise process 
may mimic pink noise due to the finite time of a run.) For 
each run, we split the 4096 noise amplitudes into overlap- 
ping bins of 256 seconds, multiplied by a Hann window, 
Fourier transformed, and squared, averaging the result- 
ing spectra in order to reduce scatter Q. 

As in the example of the figure, some runs show a 
crossover between low- and high-frequency values for a, 
which we estimated with least-squares straight-line fits 
of power at the first 100 (excluding zero frequency and 
the very lowest frequency) and last 100 frequencies. In 
many runs, low- and high-frequency alpha estimates were 
equal, within fitting errors. Table U summarizes the re- 
sults, giving in the columns labeled "ave" the means over 
all HOSE runs or all SKOV runs for the given measures. 
The differences between alphas for HOSE and SKOV, 
both low- and high-frequency, exceed several standard er- 
rors (or standard deviations of the mean, std/\/lV, where 
N is the number of runs). Moreover, the Student-i-test 
and Kolmogorov-Smirnov test show that the HOSE and 
SKOV populations differ [22|. The low- frequency expo- 
nents are more significant. The fact that these measures 
are larger for HOSE than for SKOV suggests a difference 
in long-time correlations in micromotion and is consis- 
tent with the hypothesis that non-cancerous HOSE cells 
move in a more orderly manner than cancerous SKOV. 

A non-zero aiow is indicative of long-time, "fractal" 
Q, correlation, but as Rangarajan and Ding [9] point 
out, relying on power-law behavior alone can lead to in- 
correct identification of such correlations when none ex- 
ist. Two related measures are the Hurst exponent and 
the exponent of detrended fluctuation analysis ^, HHj 
131; both methods split the time series of noise 
into bins of duration T, then determine how a measure 
scales with T. For Hurst, one subtracts the mean from 
all the data in a bin and characterizes that bin by its 
standard deviation, S. The series is integrated, and the 
minimum value subtracted from the maximum, yielding 
the range, R. For each bin, one records the ratio R/S 
and averages over bins of the same size. The procedure 
is repeated for successively larger bins (T). A straight- 
line fit to a log-log plot oi R/S against bin size T reveals 
a power law, R/S ~ , where H is the Hurst expo- 
nent. Detrended fluctuation analysis runs along similar 
lines, but within each bin one subtracts a best-flt line, 
thus detrending the data. The data in the bin are then 
characterized by standard deviation S ~ T^, where D 
is the DFA exponent. Table HTl shows the results; again, 
with high confldence (based particularly on Student's t- 
test) we can conclude that HOSE and SKOV noise come 
from different distributions. However, since the means 
are separated by less than a population standard devia- 
tion, many runs (of 4096 seconds) would be necessary to 
determine the provenance of one particular culture. 
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TABLE I; Power-spectral measures of HOSE (non-cancerous) and SKOV (cancerous) resistive and capacitive noise series. 
Shown are estimates for 1//" behavior at high and low frequencies. The means of the alphas differ by many standard errors 
(std/y/N), allowing us to distinguish the populations composed of A'^ runs, although not by enough to distinguish reliably 
a single HOSE run from a single SKOV. The _F-test and t-test give the probabilities that the variances and means of the 
distributions of values of a would differ by as much as or more than they do if the two populations had come from the same 
Gaussian distribution. KS gives the probability under the Kolmogorov-Smirnov test that the two populations' cumulative 
distributions could differ as much as they do. Small probabilities indicate that the populations differ; a probability of 0. means 
< 10"^ iV = 34 for HOSE, A = 26 for SKOV. In all cases, we apply the approximate t-test for distributions with unequal 
variances 7] . 
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0.42 


4. X 10"^ 


0.024 


capacitance 


Q^low 


0.909 
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TABLE H: Additional measures of long-time correlation in the noise time series. Hurst and detrended-fluctuation exponents. 
See Table-|T] caption for column descriptions. 







HOSE 






SKOV 


prob. 


from same distribution 




measure 


ave 


std 


std/VN 


ave 


std std/VN 


F-test 


t-test 


KS-test 


resistance 


Hurst H 


0.770 


0.0442 


0.008 


0.744 


0.0876 0.017 


3. X 10-'' 


0.17 


0.099 


DFA D 


0.854 


0.0473 


0.008 


0.806 


0.0793 0.016 


0.006 


9.8 X 10"^ 


0.057 


capacitance 


Hurst H 


0.792 


0.0474 


0.008 


0.731 


0.0886 0.017 


9. X 10"" 


3.1 X 10"^ 


0.012 


DFA D 


0.843 


0.0479 


0.008 


0.788 


0.0748 0.015 


0.017 


2.5 X 10"^ 


3.4 X 10"^ 



While aiow, H, and D were designed to estimate corre- 
lations at diverging time scales, short-time correlation is 
conveniently determined from autocorrelation, Figure[TJi, 
normalized to unity at zero lag. The lag of first zero cross- 
ing provides one natural measure of when correlation 
is lost, but since autocorrelation curves may sometimes 
reach very small, yet positive, plateaus before crossing 
zero, we also measured the lag at which the autocorre- 
lation first crosses 1/e. In a model with only short-time 
correlation, the 1/e time estimates the exponential decay 
time. However, as we discuss below, we observed signif- 
icant deviations from exponential decay, finding better 
fits to a shifted power-law decay, 

/ 1 + t \ 

autocorrelation = f — - — j . (1) 

We fit autocorrelation, for lags in the heuristic interval 
t = 1 to t = 20 seconds, using Levenberg-Marquardt 
least-squares minimization to this form to find (3q. Table 
mil summarizes results for the two crossings and Po; the 



last distinguishes the populations of HOSE and SKOV 
runs only in that (cancerous) SKOV shows much greater 
scatter in /3o, as measured by the F-test. Both crossings 
vary greatly from run to run, but the 1/e crossing in re- 
sistance and zero crossing in capacitance distinguish the 
populations of HOSE and SKOV experiments at better 
than the 95% confidence level as measured by Student's 
t-test and the Kolmogorov-Smirnov test. In particular, 
the averaged measures show shorter crossing times and 
steeper descents (/3o) for SKOV than for HOSE, again 
consistent with the hypothesis that the micromotion of 
cancerous cultures is less correlated than that of non- 
cancerous cultures. 

With the fourteen measures summarized in Tables |I]- 
IIIIl each run of 4096 seconds can be thought of as a 
point in a fourteen-dimensional space. In such problems, 
the populations might separate into two distinct, com- 
pact clusters [H, §4.2]; while the identification of clus- 
ters in high-dimensional spaces remains an open prob- 
lem in statistical research, it is common to use the 
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TABLE III: Measures of short-time correlation in the noise time series: the lag at which normalized autocorrelation (see Figure 
[ijl) falls to 1/e, the first zero-crossing of autocorrelation, and the exponent (3o from fitting the first few lags with a shifted 
power law. See the Table-|T] caption for the statistical labels. Of these measures, the 1/e crossing (in resistance) and the zero 
crossing (in capacitance) have the greatest significance in distinguishing the populations; /3o is significant only in the sense that 
the scatter is very much greater for cancerous SKOV than for non-cancerous HOSE. 
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variance-maximizing principal-component analysis intro- 
duced by Hotelling to project onto optimal subspaces, 
usually taken to be two-dimensional 16]. Figure [2] plots 
the first two principal components. While the plot shows 
a clear difference between the two populations consist- 
ing of all runs of HOSE and all runs of SKOV, overlap 
between the two clusters makes it difficult to apply the 
technique diagnostically. We found this problem to be 
generic: an exhaustive examination of pairs of princi- 
pal components (beyond the first two) produced simi- 
lar plots, with the two populations usually less distinct 
in higher-order components, while adding or subtracting 
several measures to the list of fourteen measures did not 
improve clustering. 

Thus far, the noise measures considered have shown 
that electrical noise from HOSE and SKOV experiments 
have, on average, different correlations, but they do not 
provide a reliable way to determine whether the cells 
in a single run of 4096 seconds are HOSE or SKOV. 
However, from the normalized (zero-mean, unit-variance) 
noise time series of Figure [T|d, we can extract a probabil- 
ity distribution of noise amplitudes, as in Figure [1^. Not 
surprisingly, the distribution is approximately Gaussian; 
however, subtle deviations from normal form do distin- 
guish HOSE from SKOV, even in a single run, if we ap- 
ply the Kolmogorov-Smirnov test directly to the noise. 
This test looks only at distributions of noise amplitudes, 
rather than correlations. 

To this end, we concatenate the first nine 4096-second 
HOSE resistance runs (discarding, for this purpose, the 
second halves of the 8192-second runs) to create a HOSE 
resistance reference distribution. Similarly, we create 
a SKOV resistance reference by concatenating the first 
nine 4096-second SKOV runs. Each of the remaining 
runs is tested against the two resistance reference sets. 
The same procedure is applied with capacitance data. 
In many cases, Kolmogorov-Smirnov does not show a 
match with either distribution with high probability, but 
we can compare the two probabilities: one typical HOSE 
run matches the HOSE reference with probability 0.02 



and SKOV with probability 4.7 x 10~^, so we (correctly) 
identify this run as HOSE based on the ratio of proba- 
bilities. Of 56 tested data sets (none of which went into 
the construction of the reference sets), 42 (75%) matched 
the correct reference set by this criterion, an outcome 
that would happen by chance with probability approxi- 
mately 1.2 X 10~^. We repeated the procedure using a 
second collection of four reference sets (HOSE/SKOV, 
resistance/capacitance) constructed from nine runs not 
used in making the first reference sets. Of 64 trials (none 
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FIG. 2: (Color on-line.) Projection along the first two princi- 
pal components of the fourteen-dimensional space determined 
by Tables H UlIII Blue open symbols mark the 34 HOSE runs, 
red crosses the 26 SKOV experiments. As populations, these 
two sets are distinct, but the overlap of clusters makes it dif- 
ficult to distinguish individual runs in this type of projection. 
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TABLE IV: Percentages of correct identifications. The 
Kolmogorov-Smirnov test is applied to distributions of noise 
amplitudes against HOSE and SKOV reference sets. Two 
non-overlapping choices of reference sets are used; in neither 
case did any trial run figure in a reference set against which it 
was tested. The "average" column gives percentages weighted 
by numbers of trials (16 HOSE and 12 SKOV in the first set, 
18 HOSE and 14 SKOV in the second). 





first set 
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HOSE capacitance 
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66.7 


76.5 


SKOV capacitance 


83.3 


100. 


92.3 


SKOV resistance 


66.7 


100. 


84.6 


all resistance 
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all HOSE 


75.0 


69.4 


72.1 


all SKOV 


75.0 


100. 


88.5 


all 


75.0 


82.8 


79.2 


TABLE V: Kurtosis 


averaged over 


all runs, standard devia- 


tion of kurtoses, and standard deviation of the means. 




average 


std. 


std./VN 


HOSE resistance 


74.4 


173.1 


29.7 


SKOV resistance 


3.00 


9.57 


1.9 


HOSE capacitance 


17.6 


32.2 


5.5 


SKOV capacitance 


0.94 


1.80 


0.35 



used in the new reference sets), 53 (83%) were identified 
correctly, with corresponding probability 5 x 10^*. The 
results from the two sets of trials are added and summa- 
rized in Table HVl We can reduce percentages of incorrect 
identifications by insisting on agreement between resis- 
tance and capacitance time series; this lowers the overall 
incorrect identification rate to 11.7%, with a correct rate 
of 70.0% and a "not-sure" rate of 18.3%. 

Deviations from normality in the noise-amplitude dis- 
tributions are characterized in part by kurtosis p7j . 
which is larger for HOSE than for SKOV: see TableEl A 
possible explanation is that kurtosis here is a proxy for 
correlation time: since convergence under the central- 
limit theorem is non-uniform, with a distribution ap- 
proaching a Gaussian slowly in the tails as the number 
of samples increases, a smaller population will tend to 
have a larger kurtosis than a larger one. All of our runs 
have the same number of time steps, but as we have seen, 
SKOV correlation times are shorter than HOSE correla- 
tion times; thus, a SKOV run could be said to have more 
independent time steps than a HOSE run of the same 
length and might be expected then to have a smaller kur- 
tosis. However, this cannot explain the whole effect: as 
we argue in Appendix El both kurtosis and Kolmogorov- 
Smirnov appear to be better discriminants than a direct 
measure, the 1/e crossing. This suggests that the uni- 
variate noise distribution is more than just a proxy for 
correlation time. 



IV. TWO SIMPLE MODELS 

Having motivated and interpreted our measures of 
noise in terms of short- and long-time correlations, we 
now compare our data to the simplest possible discrete- 
time models, the binary random walk with persistence 
[Tsl . displaying only short-time correlation, and a dis- 
crete fractional Brownian motion [l^ , which has corre- 
lations on all time scales. For present purposes, it suffices 
to consider only the increments rather than the walks 
themselves; that is, we compare to Figure [TJd, not Figure 
[Da. 

First, consider the increments of a discrete random 
walk with persistence. Let the increment at time jAt, 
where At is the time step, be Xj, drawn from { + 1, — 1}- 
Then Xj+i — Xj with probability a and Xj^i = —Xj with 
probability 1 — a; one recovers the usual discrete binary 
random walk for a = 1/2. Since we think of this pro- 
cess as approximating a continuous one, and there is no 
natural way to take the limit At — > for anticorrelated 
increments, we restrict 1/2 < a < 1. For convenience, we 
set At = 1. A simple inductive argument shows that 

(xoxn) = {2a - 1)" = exp(-n/T) , (2) 

where the correlation time r = — l/ln(2a— 1). For times 
much larger than r, this Markov process looks like an or- 
dinary binary random walk with a rescaled time, and by 
the usual arguments [13] , the power spectrum approaches 
white noise, i.e., it is independent of frequency. However, 
for a finite run, the power spectrum may mimic corre- 
lated (pink) noise even, surprisingly, for a t as short as 
4 in a run as long as 4096, as in Figure [5^. However, 
the random noise levels off noticeably at low frequencies, 
while the experimental data (Figure [3)3) appear to follow 
a 1//" power law to the lowest frequencies [23|]. This 
supports the presence of correlations at all time scales. 
The shortness of the low-frequency plateau in Figure[3^ is 
misleading. To see more of the flat part of the spectrum, 
finer frequency resolution is necessary. Taking larger win- 
dows, we can (at least for a run longer than 4096) extend 
the graph many decades to the left and verify that the 
spectrum remains flat (white), but at the cost of greater 
scatter. Parts (c) and (d) of the figure show autocorre- 
lation for random noise and experimental data with fits 
to exponential decay (dotted) and the shifted power law 
(H]) (solid). The two fits fall on top of one another for 
the process satisfying That exponential decay does 
not approximate the experimental data as well as the 
power law corroborates the hypothesis of longer-than- 
short-time correlations. 

Mandelbrot and van Ness introduce the notion 
of fractional Brownian motion with correlations between 
increments separated by arbitrary time differences and 
with a 1//" power spectrum. Rangarajan and Ding 
describe a particularly simple way of generating a 
time series of increments with such properties: start 
with a Gaussian-distributed uncorrelated time series 



6 




persistence times ranging from 2 at the left to 7 at the 
right. Measured values of a follow the same prescription 
as used above, although as noted earlier (Figure [3|) , the 
fits fail for low frequencies; indeed, every a should be 
zero. Hurst estimates range from 0.45 to 0.67; the true 
value in every case should be 1/2. As discussed by Ran- 
garajan and Ding, the discrepancies between measured 
Hurst and Hurst estimated from measured a are large. 
In the middle and at the bottom are plotted our exper- 
imental data (HOSE o, SKOV x). Agreement between 
the exponents H and a is generally not as good as for 
the long-range-correlated processes but not so poor as for 
the short-time-correlated random walk. On average, the 
experimental points lie closer to the former than to the 
latter. We interpret this result as supporting the exis- 
tence of correlations on, at the very least, many different 
time scales. A model of cell motion will need to explain 
both the short-time and long-time correlations we have 
observed. 



FIG. 3: The increments of a finite random walk with per- 
sistence (left) may mimic certain aspects of the experimental 
data (right), but with notable differences. The random pro- 
cess has a — 0.8894, so an exponential decay time r = 4.00. 
The experiment is a typical capacitance noise time series of 
HOSE, with a measured 1/e crossing of 5.7. (a) and (b) show 
the best-fit lines to the first 100 points (excluding zero and 
the lowest frequency) of the power spectrum; both give slopes 
~ —1.0, but the random data level off noticeably at low fre- 
quencies, as would be expected of white noise. Autocorrela- 
tion curves (c) and (d) show fits to exponential (dotted line) 
and shifted power-law ^ (solid) decays. For the random 
noise, the two fits fall on top of one another, but for the 
experimental data, a power law fits better than exponential 
decay. 



{xj}, Fourier-transform, multiply by /~"/^, and Fourier- 
transform back. The resulting process has a Hurst expo- 
nent given by 



H ^ (1 +a)/2 



(3) 



Determination of exponents a and H is subject to the 
usual numerical vicissitudes, but Rangarajan and Ding 
argue that true long-ranged processes should satisfy ^ 
at least approximately. 

Figure m plots fractional discrepancies between ([3]) and 
measured Hurst exponents as functions of measured spec- 
tral exponents a. At the bottom are plotted artificially- 
generated long-time-correlated data following the pre- 
scription of Rangarajan and Ding (plotting symbols -I-); 
the measured exponents a are always close to the known 
values, so the measurement errors occur in estimat- 
ing H. We note a systematic trend toward larger er- 
rors away from a « 0.5, but generally the errors stay 
small. At the top of the graph (plotting symbols o) 
are artificially-generated random walk increments with 
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FIG. 4: (Color on-line.) Fractional discrepancies between 
i^prcd given by (|3} and measured Hurst exponent as func- 
tions of measured spectral exponent a. Near the bottom, 
plotted with large -I- symbols, are artificially-generated data 
with known long-time correlations. At top are generated data 
(large o) from random- walk increments with persistence times 
ranging from 2 (smaller values of a) to 7 (larger values). In 
the middle are experimental results for HOSE (blue o) and 
SKOV (red x). Most of the experimental data look more like 
the correlated data than the uncorrelated, but a few overlap 
with uncorrelated noise; all of these are SKOV. 
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V. APPLICATIONS 

We have demonstrated that electrical-noise measure- 
ments on human ovarian surface epithelial cells can dis- 
tinguish cancerous and non-cancerous cultures. This is 
not intended as a diagnostic tool; for one thing, it is easier 
to distinguish them under a microscope. We find it is also 
possible to distinguish HOSE from SKOV based purely 
on average electrical resistance or capacitance. Our main 
focus has rather been on developing statistical tools with 
which to test more sophisticated statistical-mechanical 
models and in developing a database of characteristics 
of many different cell types, for which a single measure- 
ment {e.g., average electrical resistance) will surely be 
inadequate. 

One characteristic of malignant cancerous cells is their 
ability to invade tissue in disregard of clues from their 
neighbors 21|; our observation of shorter correlation 
times in cancerous cultures is consistent with the pic- 
ture of cancerous cells moving in a less regulated man- 
ner. Now that it has been established that different cell 
types generate distinguishable noise patterns, future re- 
search in this area will focus on the development of realis- 
tic models of cellular motility for healthy and malignant 
cells. 
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APPENDIX A: COMPARING DISCRIMINATORS 

We claim no originality to the following elementary 
application of statistics but could not find a textbook 
discussion of quite this point. Given two distributions, 
A and B (for instance, the kurtoses of HOSE data sets 
and those of SKOV), assumed to be Gaussian and char- 
acterized by means < and standard deviations 
cta, ctb, there are several choices of where to place a di- 
viding point a;o so as to identify all x < xq as belonging 



to population 1 and all a; > a;o to population 2. One 
natural choice is to pick xq so that the expected rates of 
correct identification of the two populations will be the 
same, i.e., that xq — should be the same multiple of 
a A as ^iB — xq is of ctb, or 



Xo 



<7A + CTb 



(Al) 



Any other choice will decrease the expected rate of in- 
correct identification of one population at the cost of in- 
creasing the other. A second plausible choice is to seek to 
maximize the sum of the expected correct identification 
rates, 



1 1 r 



Xo - flA 



CTA 



^/2 



11 r f f^B - Xq 



(A2) 



Cr = - + -erf 



2 2V <7bV2 



it is easy to show that the separatrix xq is then 



2^0 = ■ 



fJ.B<y\~crB (^fJ:A<yB±<yA\/ (AiA-MB)^+2(cr^-(T|) ln((TA/o-B) ) 

^ 5 

'^A-'^B 

(A3) 

(One root maximizes C^i + Cb- Note that (jA3| reduces 
to [fiA + fJ'B)/'2 when cta = as.) A third natural choice, 
maximizing the product Ca C b , requires numerical solu- 
tion. Of course, a more complicated risk function could 
apply, for instance in medical diagnosis, where a false 
negative is much worse than a false positive. 

To compare the predictive values of three of the statis- 
tical measures developed in the text, 1/e crossing from 
Table IHIl kurtosis from Table |Vl and the Kolmogorov- 
Smirnov test of Table IIV| we apply the simplest sepa- 
ratrix, (|A1[) to the means and standard deviations esti- 
mated for the first two. (This choice is motivated by the 
similar correct-identification percentages for HOSE and 
SKOV in Table [TVl but as an alternative to (|A2p . using 
the actual data sets gives comparable answers). Then 
the expected correct-identification rate (|A2[) for 1/e as 
a discriminant is 62% and that for kurtosis 67%. These 
rates are both lower than the 79% (Table IIV|) for the 
Kolmogorov-Smirnov test applied to the noise distribu- 
tion, undermining the idea that the deviation of this dis- 
tribution from normal form is strictly a proxy for corre- 
lation time. 
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